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Abstract 

We use the line digraph construction to associate an orthogonal matrix with 
each graph. From this orthogonal matrix, we derive two further matrices. The 
spectrum of each of these three matrices is considered as a graph invariant. For the 
first two cases, we compute the spectrum explicitly and show that it is determined 
by the spectrum of the adjacency matrix of the original graph. We then show 
by computation that the isomorphism classes of many known families of strongly 
regular graphs (up to 64 vertices) are characterized by the spectrum of this matrix. 
We conjecture that this is always the case for strongly regular graphs and we show 
that the conjecture is not valid for general graphs. We verify that the smallest 
regular graphs which are not distinguished with our method are on 14 vertices. 
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1 Introduction 



Graphs are often conveniently represented using matrices, for example, the adjacency 
matrix, the Laplacian matrix, etc. [5 a . Many important properties of a graph are encoded 
in the eigenvalues of the matrix representation. However, eigenvalues generally fail to 
separate isomorphism classes. In this paper, we consider some matrix representations 
inspired by the notion of coined quantum walks Since strongly regular graphs give 
rise to relatively large sets of non-isomorphic graphs which are cospectral with respect 
to the commonly used matrix representations, we use these graphs as a testing ground 
for the representations that we define. Let G be a graph and let U be the orthogonal 
matrix inducing a coined quantum walk on G, where the coins are Grover matrices. For 
every given k, we define a digraph, D, with an arc if and only if U^- > 0. For 

k = 1, we express the adjacency eigenvalues of D in terms of the ones of G. When 
considering strongly regular graphs, we describe how to construct directly D from G. For 
strongly regular graphs, we conjecture that the eigenvalues of the adjacency matrix of D 
distinguish G from its cospectral mates for k = 3. We verify this conjecture numerically 
for all strongly regular graphs up to 64 vertices contained in the tables of Spence ^1] 
(the graphs srg(16, 9, 4, 6) were obtained from the tables of McKay ;9j). By providing 
counterexamples, we show that our method fails to distinguish general graphs. Hopefully 
these counterexamples will help understanding what graphs are cospectral with respect 
to the representations described here. 

The remainder of the paper is organized as follows. In Section 2, we recall the defini- 
tion of the orthogonal matrix U and give a formula for its eigenvalues. We define digraphs 
derived from powers of U, and propose using these digraphs to distinguish the original 
graph from its cospectral mates. In Section 3, we focus on strongly regular graphs. We 
describe some structural properties of the digraph D obtained from U 3 . In Section 4, we 
list the sets of strongly regular graphs for which the eigenvalues of D successfully distin- 
guish a graph from its cospectral mates. We conclude with counterexamples involving 
general graphs. 

Recently, dynamical processes based on quantum evolution have been considered for 
attacking the graph isomorphism problem jHJEl- in this context, Shiau et al. jT^] gave 
evidence that quantum walks alone are not successful. The setting they describe is as 
follows. Let {\j) : 1 < j < n} be an orthonormal basis of a Hilbert space is a ray 
and (i\ is the linear functional which maps each \j) to the usual inner product, which is 
denoted by Given a graph G, let us define the Hamiltonian H — — ^ M{G) a j,c\cb, 

where the operator c\cb is given by {i\c\cb\j) = 5i, a 8b,j and M(G) is the adjacency matrix 
of G. The evolution for a time t of the initial states {\ipj(0)) = \j) : 1 < j < n} is 
governed by the Schrodinger equation i d ^ 3 Jf^ = H\ijjj(0)}. In |13j . it was shown that 
the matrix whose entries are Ojj = (ipi(0)\ipj(t)) do not help in distinguishing pairs of 
strongly regular graphs with the same set of parameters. However, by making use of some 
techniques described in (TT). it was also given an alternative method which distinguished 
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numerically all known strongly regular graphs up to 29 vertices. It would be interesting 
to investigate potential connections between [T3] and the present paper. 



2 Representing graphs with orthogonal matrices 

Let G = (V, E) be a simple undirected graph (that is G is loopless and without multiple 
edges). Let D G = (V, A) be the directed graph obtained from G by replacing every edge 
G E(G) with the pair of arcs and (j, i). The degree of a vertex i is denoted by 
d(i). A graph is said to be fc-regular if every of its vertices has degree k. Unless otherwise 
stated, in this paper we will consider only graphs with minimum degree 3. 

Definition 2.1 Given a graph G, we denote by U(G) the matrix defined as follows: 



The matrix U(G) induces a coined quantum walk on G, where the coins are Grover 
matrices [I]. Here are some basic observations about U(G): 

• If \E(G) | = m then U(G) is 2m x 2m and the vertices of degree zero do not contribute 
to the dimension of Uq- 

• Given i G V(G), if d(i) = 1 and {i, j} G E(G), then U(G)ufi = 1, and this is 
the unique non-zero entry in the row indexed by (j, i) and in the column indexed 
by 

• Given j G V(G), if d(j) = 2 and {i,j},{j,k} G E(G), then U{G) {i>j)m = 
u ( G )(k,j),(j,k) = and U{G)^ j )^ hk) = U(G)( k ,j),V,i) = 1- So, if G is 2-regular, 
the matrix U(G) has exactly a one in every row and every column and it is then a 
permutation matrix. 

• Since G is undirected, ^"(G0(ij),(j,i) 7^ and U{G){j y i)^j) ^ for every G 
E(G). However, we have assumed that the minimum degree of G is 3, and then 
U(G) is not symmetric. In fact, for three distinct vertices G V(G) such that 
{M'}> {j 5 e we have U{G)^) t (j t i) ^ 0, but U(G)<j tV) ^ = since z 7^ /. 

The support of an n x n matrix M, denoted by M, is the n x n (0, l)-matrix with 
elements defined as follows: 



Let M(G) be the adjacency matrix of a graph (digraph) C The line digraph of a 
digraph D, denoted by LD, is the digraph defined as follows: V(LD) = A(D) and 



for all (z,j),(k,l)eA(D G ), 






otherwise. 
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(k, I)) G A(L D) if and only if j = k. If d(i) > 2 for every i G V(G) then one can 
verify that U(G) = M(LD G ) [12]. 

Graphs (or digraphs) G and H are isomorphic (written G = H) if there is a permuta- 
tion matrix P such that M(G) = PM(H)P~ 1 , where M(G) and M(H) are the adjacency 
matrices of the graphs G and H, respectively. Harary and Norman [7j proved that if 
D and F are digraphs without sources or sinks then L D = L F if and only if D = F 
{i G V(D) is a source if there is no j such that (j, i) G A(D); a smA; if there is no j such 
that ^ A(Z})). Since .D^ is without sources or sinks, it follows that G = H if and 
only if U(G) = QU(H)Q~ 1 , for some permutation matrix Q. According to Definition 1, 
this fact is sufficient to ensure that G = H if and only if U(G) = QU(H)Q~ 1 , for some 
permutation matrix Q. In words, U(G) "faithfully" represents G. 

The spectrum of a matrix M is the multiset of the eigenvalues of M, and is denoted 
by Sp(M) = {[Aj] mi : 1 < i < m}, where Ai, A 2 , are the eigenvalues of M and 
m!,m 2 , ■■■,rrik their respective multiplicities. Let T{G) be the matrix with r/'-th entry 



We can obtain Sp(£7(G)) from Sp(T(G)), as shown in the next result. 

Proposition 2.2 Let G be a graph on n vertices and m edges. The matrix U(G) has 2n 
eigenvalues of the form 



where At is an eigenvalue of the matrix T(G). The remaining 2(m — n) eigenvalues of 
U(G) are ±1 with equal multiplicities. 



Proof. Let M = M(G), U = (G) and T = T(G). Since G is undirected, we have 
Mf - = M ifj M jti = M itj . By Definition 1, 



Let u be an eigenvector of U with eigenvalue A (notice that the entries of u have two 
indices as it is for the rows and the columns of U). By making use of this equation, we 
can write 







k,l 
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Let t be an eigenvector of T with eigenvalue Ay. Define 



u 



(id) 



d(j) 



d(i) ' 



Then 



2^1 
d(j) 



i) 



Mi d 



i,J 3 



(2 A T - A*) - M i; 



U 



Since U is real-orthogonal, its eigenvalues have magnitude 1. It follows that 



An 



(iJ) 



A 



Midi. 
— ^(2A*A 



(A 



*\2\ 



A*M 



Therefore, u is an eigenvector of U if 2A*At — (A*) 2 = 1 or A + A* = 2At (A* is the complex 
conjugate of A). It is easy to verify that the remaining eigenvalues are ±1. □ 

Two graphs G and H are said to be cospectral (or isospectral), with respect to a given 
matrix representation, if the spectra of the matrices representing the graphs are identical 
(see, e.g., JTOJ ) - For example, let G and H such that 



M(G) 



and 



M(H) 
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Then Sp(M(G)) 
of Proposition 



Sp(M(if)). The following corollary is a straightforward consequence 

Corollary 2.3 Given graphs G and H , Sp(C/(G)) = Sp(U(H)) iff Sp(T(G)) = Sp(T(#)). 

If G = H then Sp(M(G)) = Sp(M(if)), but the converse is not necessarily true. This 
is the case for the graphs G and H considered in the above example, since Sp(M(G)) = 
Sp(M(#)) = {[-2] 1 , [0] 3 , [2] 2 } but, of course, G£H. Also, Sp(T(G)) = Sp(T(£T)) does 
not imply that G = H. For instance, if G is a fc-regular graph then Sp(T(G)) = {[^Aj] mi : 
1 < % < k}, where Sp(M(G)) = {[Aj] mi : 1 < i < k}. So, for two fc-regular graphs G 
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and H, Sp(M(G)) = Sp(M(#)) if and only if Sp(T(G)) = Sp(T(#)), or, equivalently, 
Sp(C/(G)) = Sp(U(H)). Two non-isomorphic /c-regular cospectral graphs G and /J are 
then not distinguished by the spectra of U(G) and U(H). 

Now, let M and N be matrices such that Sp(M) = Sp(iV) = {[A;] m * : 1 < i < n}. For 
any natural number p, Sp(M p ) = Sp(iV p ) = {[Af]™ 1 : 1 < % < n}, but it is not necessarily 
the case that Sp( M p ) = Sp(iV p _). This is again the case of the graphs G and H in the 
above example: 



M(Gf 
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M{Hf = 



Sp( M(G) 2 ) = {[0] 3 , [2] 2 } and Sp( M(ff) 2 ) = {[0] 3 , [1], [4]}. The entry M(G) P 5j equals 
the number of paths of length p between vertices i and j. Note that M{G)\ 1 = 2 and 
1. This implies that matrix M(G) P contains generally less information about 



M(G)< 

the structure of a graph G, since M(G) P . = 1 if and only if the number of paths of 

length p between vertices i and j is non-zero. Despite this fact, it is legitimate to ask 
the following question: given graphs G and H such that Sp(M(G)) = Sp(M(iJ)) what 
conditions imply that Sp(M(G) p ) ^ Sp(M(H) p ), for some pi If a graph G with diameter 
p is connected and non-bipartite then M(G) P = J, where J denotes the all-ones matrix. 
As a consequence, if two graph G and H have diameter p then Sp(M(G) q ) = Sp(M(H) q ) 
for every q > p, even if Sp(M(G)) ^ Sp(M(#)). The spectra Sp( M(G) 9 ) and Sp(M(ff)g) 
are indeed of no use in distinguishing G and if if these graphs have diameter 2. 

Given a graph G, since the matrix U (G) is real-orthogonal, the support of U (G) p , for 
some p, is presumably different from the support of M(L Dg) p ■ This is because of the 
contribution of the negative entries in U(G). On the basis of this observation, we define 
the following matrix. 

Definition 2.4 Given a graph G, we denote by S + (U(G) P ) the matrix defined as follows: 

Suppose that Sp(M(G)) = Sp(M(#)) and G ^ H. What conditions need G and 
# to satisfy in order that Sp(S + (U(G) p )) ^ Sp(S + (U(H) p )) for some p > 2? Strongly 
regular graphs seem to provide a good testing ground to approach this question. Strongly 
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regular graphs have been investigated in many different contexts, including group theory, 
algebraic graph theory, design of experiments, finite geometries, error-correcting codes, 
etc. (see, e.g., |21E1)- The fastest known algorithm for testing the isomorphism of strongly 
regular graphs on n vertices was designed by Spielman [T3] and runs in n 0(nl 3 logn ^ number 
of steps. A strongly regular graph with parameters (n,d,r,s) (for short, a srg(n, d, r, s)) 
is a (i-regular graph on n vertices such that any two adjacent vertices have exactly r 
common neighbours and any two nonadjacent vertices have exactly s common neighbours 
Strongly regular graphs are interesting in our context principally for the following 
two reasons: 

• For the strongly regular graphs G and H with identical sets of parameters, we 
have Sp(M(G)) = Sp(M(if)) even if G ¥ H. In particular, the adjacency matrix 
of a srg(n, d, r, s) has exactly three eigenvalues, [d] 1 > [e + ] m+ > [e~] m , such that 
e + = |(r— s+\/A) and e~ = \{r— s— \/A), with {m + , m~} = \{n— 1 ± 2k +(n-^s-r) ^ 

and A = (s - rf + 4 (d - s). 

• The diameter of a connected strongly regular graph is 2 (see, e.g., [3]). 

3 On the spectra of S + {U{G)), S + (U(G) 2 ) and S + (U(G) 3 ) 

In this section, we give formulas to compute the eigenvalues of S + (U(G)) and S + (U(G) 2 ) 
from the eigenvalues of M(G), when G is strongly regular. It follows that neither 
S + (U (G)) nor S + (U(G) 2 ) are of any help in distinguishing G. The eigenvalues of S + (U(G) 3 ) 
do not depend entirely on the eigenvalues of M(G) and it seems difficult to give a for- 
mula for Sp(S + (U(G) 3 )) in terms of properties of G. On the other hand we show how to 
construct S + (U(G) 3 ) from G. We conjecture that Sp(S + (U(G) 3 )) distinguishes G from 
its cospectral mates. 

Given a digraph D, let P(D,x) = det(xJ — M (D)) be the characteristic polynomial 
of M (D). It is known that |Hl E3 

P(LD,x)=x\ A ^\-\ y ^P(D,x). 

Then 

Sp( U(G) ) = Sp(M(G)) U {[o]nk-(i+m++ m -)^ 

The following two propositions show that the spectra of S + (U(G)) and S + (U(G) 2 ) 
are determined by the parameters of G. Therefore the spectra of these matrices do not 
characterize G itself. 

Proposition 3.1 Let G be a k-regular graph on n vertices. Let e be an eigenvector of 
M{G) with eigenvalue Am and let u be an eigenvector of S + {U{G)) with eigenvalue X. 
Then S + (U(G)) has 2n eigenvalues of the form 
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with eigenvectors having entries 

",, = M(G) id ej - {(:),.,<,. 

The remaining n(k — 2) eigenvalues of S + (U(G)) take the values ±1. 
Proof. Observe that 



S + (U(G\ iJUk>l) ) = {l-5 hl )M h3 M Kl 5 hk . 



We then have 



(a,b),(x,y) ) u x,y 
= M a)b ^ M b,yUb,y ~ M a , b U bja 

y 

= M afi M b, y e y - jM a , b e b 
y 

M b,y - M atb e a + jM a , b e b 

y 

= M a , b e b (\ - (k - 1)/A) - M a>b e a 
= A ^M a ^ b e b - jM a , b e a 

Then these eigenvalues and eigenvectors account for 2n of the eigenvalues and eigenvectors 
of S + (U(G)). It is easy to verify that the remaining n(k — 2) eigenvalues take the values 
±1. □ 

Proposition 3.2 Let G be a k-regular graph on n vertices. Let e be an eigenvector of 
M(G) with eigenvalue \m and let u be an eigenvector of S + (U(G) 2 ) with eigenvalue A. 
Then S + (U(G) 2 ) has 2n eigenvalues of the form 



A = -f- + 2 - k ± i\ M sjk-l- A^/ 4 
with eigenvectors having entries 



U i,j = M \^ r )i,j e j - Y^{k~^V) j6i ' 

The remaining n{k — 2) eigenvalues of S + {U(G) 2 ) take the value 2. 
Proof. Observe that 

S + {U{G) 2 ) = M hJ M hk M k ^25 hk 5 Jtl + 1 - S i>k - 5 jtl ). 
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We then have 



We can also write 



If 

then 



S + (U(G)^ b) ^ y) )u X: y 

= ^2 M a:b M b:X M Xt y(25 atX 5 bi y + 1 - S a>x - 5 by y 

A M (A*-2 + fc) ^ 
tvl x,y e y k — 1 m x,y e x 

M a , b e b {\ 2 M - A*) - M a , b e a 

(2-fc)(A* -2 + fc) + A^(fc-l) 

X M (k-l) 
M a , b e b \ - M a , b e a 
(A* - 2 + fc) (2 - fc + A - 2 + fc) 
A M (fc-l) 



^ 5 + (f/(G')( ai6 ) i ( a;i j / ))M( ;E) y) 

(z,2/) 

y u x,y 

x,y 

- M a , b J2 M b:a M a:y u a , y - M af) J2 M b:X M xfi u x , b . 
y % 

T, v M a, y u a ,y = and V M xfi u x , b = 



(x,y) 

and u is an eigenvector with eigenvalue 2, provided that M a b = and then u a;b = 0. □ 

The following observation allows us to construct S + (U(G) 3 ) directly from a strongly 
regular graph G, without the need of first constructing U (G) . 

Proposition 3.3 Let G be a srg(n, k, r, s) . Then S + (U (G) 3 { ^ ^ m )) = 1 if and only if one 
of the following conditions holds: 

1. i = m, j 7^ I and 

k 2 

s + (r - s) M(G) hl - k + — > 
(which always holds if i — m, j ^ I and k > A); 
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2. i = l,m^j and M(G) hm < f ; 

3. i = I and m = j; 

I i^i >m = j an d M(G) it i < f ; 
5. i ^ I, i ^ m, j ^ I, j ^ m and 

2 [s + (r _ S )M(G) W ] > M(G) ii i + M{G) hm . 

Proof. Let G be a srg(n, A;, r, s). Let A = M(G) and U = U (G). We define the amplitude 
of a given path to be the product of the entries of U corresponding to the transitions along 
the path. The entry S + (U (G) 3 ^ ~ „ ,) is given by the sum of the amplitudes of all possible 
paths of length three from (i,j) to (1,171). In what follows we count all these paths to give 
conditions for S + (U(G) 3 i ^ ^ m ^) = 1. We begin by observing that i ^ j and I ^ m since 
G is a simple graph. The cases in the following table completely describe the possible 
ways in which the vertices /, and m can be identified with one another: 



i = m(^ I) 


i = l(^ m) 


i ^ m and % ^ I 


Case A: j = I 
Case B: j ^ I 


Case C: j = m 
Case D: j ^ m 


Case E: j = m(^ I) 
Case F: j — l(^ m) 
Case G: j ^ m and j ^ I 



For each case we count the possible paths below and thus determine y m y Note 
that all paths of length three must be of the form 

(hj) -»• (j» -»■ -> (l,m), 

where x G {i, m} or x is some vertex other than i,j, I or m. 

Case A (i — m and j ^ I) There is the unique path that takes place only on the 
vertices I, m}: 

(hj) -»• (j,i) = (j,m) -> (m,l) -> (l,m), 
and has amplitude |(| — l) 2 . The other paths are of the form 

(hj) -»■ (j» -»■ OM) -> (i,»ri), x i {i,j,l,m}, 

which have amplitude p-. If {j, 1} G E(G) there are r — 1 such paths, if {j,l} £ E(G) 
there are s — 1. Thus, 

= | { (f - x ) 2 + ^ + <r - ^ - !]} ■ 

Case B (i = m and j = /) There is the unique path that takes place only on the 
vertices {i,j, I, m}: 

(hj) -> (j,i) -»• (i, j) -> CM) = 
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which has amplitude (| — l) 3 . There are k other paths of the form 

(hj) -> U,x) -> -> = (Z,m), x £ {i,j,l,m}, 
which have amplitude p. Thus, 



2 

'(i,j),(l,m) — 



Case C (« = / and j ^ m) There are no paths which take place only on the vertices 
{i,j,l,m}. There are paths 

-»• (J» -»■ -> (Z,m), x £ {i,j',Z,m} 

of amplitude p. If {j, m} e -E'(G') there are r — 1 such paths, if {j,m} ^ E(G) there are 
r. Thus, 

U {i,j),(l,m) ~ p I y A 3,rn 

Case D (i — I and j = m) There are no paths which take place only on the vertices 
{i,j, l,m}. There are r paths 

(hj) -»• (j» -»• -> (Z,m), x ^ {i,j,/,m} 
with amplitude p. Thus, 

3 8 
^(ij),(i,m) = 

Case E (i ^ m,i ^ I and j = to) If {«,/} ^ E(G) then there are no paths which 
take place only on the vertices l,m}. However, there are r paths 

(hj) -»■ (j» -»■ (z,Z) -> (l,m), x £ {i,j,l,m} 
with amplitude p. If {«, /} G -E'(G') there is the unique path 

(hj) -»• -»• (M) -»• 
with amplitude p(| — 1). When {i, 1} e -E'(G') there are also r — 1 paths 

(m) -»■ (j» -»• (z,Z) -»• (l,m), x £ {i,j,l,m} 
with amplitude p. Thus, 

L, (tJ),(i,m) - p ^ ~ A ^ 

Case F (i ^ m,i ^ I and j = /) There are two paths 

(hj) -> CM) -> (i, j) -> (j,™) = (^ m ) 
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and 

each of amplitude |(f — l) 2 - There are also k — 2 paths 
with amplitude — 1). Thus, 

U (i,j),(l,m) = 0- 

Case G (i ^ m,i ^ l,j ^ I and j 7^ Z) We consider the cases {j, /} G 
and {j,l} £ E{G) separately. Firstly, we consider the case {j, 1} G E{G). There are 
r - Ai t i - Aj^ m paths 

(hj) -> (j> x ) -> -> (Z,m), a; ^ {i,j,Z,m} 

of amplitude p. There are also A^j + A,- )m paths of amplitude p(| — 1). If they exist, 
these are 

and 

CM) (*\0 

Now, if {j, /} ^ -E(G) the analysis is identical but with s replacing r. Thus, 



Notice that the conditions in the statement of the proposition are mutually exclusive. 

■3 

(i,j),(l,m) 



Each case corresponds to Uf- ,-, (l ) > for one of the cases A-F above. The condition 1 



corresponds to the case A, 2 to C, 3 to D, 4 to E and 5 to F. We have U 3 { ^ < for 
all the other cases. □ 

We conclude with the following conjecture. 

Conjecture 3.4 Let G and H be strongly regular graphs with the same set of parameters. 
Then G = H if and only if Sp(S + (U (G) 3 )) = Sp(S + (U(H) 3 )). 

If the conjecture is true then isomorphism of strongly regular graphs can be tested in 
polynomial time. 



4 Computations 

The table below contains the parameters of the strongly regular graphs for which we 
have tested the validity of the conjecture. We have tested also the complements of these 
graphs. In all cases the conjecture resulted to be valid. The number of non-isomorphic 
srg(n, k, r, s) is denoted by N(n,k,r,s). The graphs srg(16, 9, 4, 6) where obtained from 
all other graphs where obtained from [T4] : 
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(n, k, r, s) 


AT / 1 \ 

J\{n, k, r, s) 


[n, k, r, s) 


A T i 1 \ 

J\{n, k, r, s) 


(16,6,2,2) 


2 


(35,18,9,9) 


227 


(16,9,4,6) 


2 


(36,14,4,6) 


180 


(25,12,5,6) 


15 


(36,15,6,6) 


32,548 


(26,10,3,4) 


10 


(40,12,2,4) 


28 


(28,12,6,4) 


4 


(45,12,3,3) 


78 


(29,14,6,7) 


41 


(64,18,2,6) 


167 



The method that we have described appears then to be successful for strongly regular 
graphs. In general, there are cases in which the method fails. We have verified that 
Sp(S ,+ (f/(G) 3 )) distinguishes all regular graphs up to 13 vertices. However, we have 
found 4-regular graphs on 14 vertices which are not distinguished by our method. Two 
graphs for which the method fails are drawn in the figure below: 




The computer programmes used for our calculations can be downloaded from 
www-users . cs . york . ac . uk/~wilson/qwalks . html 
or asked directly to the authors. These are MATLAB and C codes. 
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